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Imposing twisted boundary conditions on the fermionic fields is a procedure extensively 
used when evaluating, for example, form factors on the lattice. Twisting is usually performed 
for one flavour and only in the valence, and this causes a breaking of unitarity. In this work we 
explore the possibility of restoring unitarity through the reweighting method. We hrst study some 
properties of the approach at tree level and then we stochastically evaluate ratios of fermionic 
determinants for different boundary conditions in order to include them in the gauge averages, 
avoiding in this way the expensive generation of new conhgurations for each choice of the 
twisting angle, 0. As expected the effect of reweighting is negligible in the case of large volumes 
but it is important when the volumes are small and the twisting angles are large. In particular we 
hnd a measurable effect for the plaquette and the pion correlation function in the case of 0 = 7r/2 
in a volume 16 x 8^, and we observe a systematic upward shift in the pion dispersion relation. 
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1. Introduction 

Non-Periodic Boundary Conditions (NPBCs) for the fermions are typically employed on the 
lattice in order to obtain a fine resolution of momenta in the spatial directions. Twisting [1] amounts 
to imposing 




, ju = 1,2,3 
y/(x) , IX =0 


where Nji is the lattice extent in direction fl and [0,27r] is an angle. Alternatively one can 
introduce a constant U( 1) interaction with vanishing electric, magnetic field and elecfric pofenfial 
buf consfanf vecfor pofenfial [2] (see also [3] for fhe equivalence of fhe fwo procedures). The above 
exfra inferacfion is implemenfed by Iransforming fhe sfandard QCD links (U) in fhe following way 


\Uo{x), Ai=0 


( 1 . 1 ) 


The modificafion of fhe boundary condifions proved fo be beneficial for: 

• Form factors: one can scan values of fhe exchanged momenfa in fhe scaflering process wifh 
a very fine resolution, in order fo defermine more accurafely fhe form factors on fhe laffice. 
An example is fhe semi-lepfonic decay Kij [4], used to exfracf fhe CKM elemenf Vus- 

• Mafching befween HQET and QCD: in fhe Schrddinger funclional, one compufes finile vol¬ 
ume observables for differenl values of d in fhe valence [5]. In parficular in such a sefup 
reweighfing fo a unifary formulalion is expecfed to be efficienl as fhe volumes considered are 
rafher small (fhe reweighfing facfors are exfensive quantifies). 

• Dispersion relafion: which is used here more as a consisfency check, as done also in [1]. 

Usually one performs fhe fwisfing only in fhe valence, which causes breaking of unifarify. This 
is simply undersfood since fhe procedure yields differenl propagafors for fermions in fhe valence 
and in fhe sea. One way to overcome fhis problem is fo perform direcf simulations wifh fermionic 
NPBCs for each value of d. Thai Ihough, would clearly be too expensive from fhe compulafional 
poinl of view. However, Ibis kind of breaking of unifarify is expecfed fo be a finile volume effecl, 
in ;^-PT for example lhal is fhe case [3]. This suggesfs lhal where fhe effecl is large, reweighfing 
may provide a reliable way to restore unifarify. 


2. Reweighting 

Lei us suppose we wanl fo conned simulalion resulls for a choice of bare paramelers a = 
{15,mi,m 2 ,... ,mnpdii ,...} to a (slighlly) differenl sel b = {P',m[,m 2 ,.. .,m{^,d '^,...} of param¬ 
elers. This can be achieved by numerically compuling on fhe a-ensemble fhe reweighfing faclor 
^ab = Pb/Pa, which is fhe ralio of fhe fwo probabilily dislribufions and if is clearly an extensive 
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quantity, Pa[U] = n”£i det(D[i7,6] +m,), where we have explicitly indicated the depen¬ 

dence of the Dirac operator on the twisting angle. In this way the expectation values on the b- 
ensemble are calculated as 


{^)b 


{ffWab)a 

{Wab)a 


with ^ the observable after Wick contractions, and the expectation value over the set of bare 
parameters a. 

By choosing to change only the boundary conditions from one bare set to the other we arrive at the 
following expression of the reweighting factor 


We=det{Dw[U,G]D„^[U,0])=det{Dw[^,0pwP,0]), 


where Dw is the massive Wilson operator. We therfore need a stochastic method to estimate (ratio 
of) determinants. For this purpose we use the following integral representation of a normal matrix 
determinant with spectrum A (A) that holds if and only if the real part of each eigenvalue is larger 
than zero [6] 

= y^^[T7]exp (-t]'*‘At 7) < oo MeA (A) > 0. (2.1) 

If it is so then the stochastic estimae converges. We use as probability distribution p(t]) of the rj 
vectors a gaussian one, then the determinant reads 


1 

detA 




It should be noted that if we require, for the case of an hermitian matrix, the existence of all gaussian 
moments then all the eigenvalues must be strictly larger than one 


/„-2ri'’Ari\ r i 

“ J < “ ^ ^ (•^) > 2' 

= J si[n]exp[-n^NA-{N-i)i]n] <» <=> a ( a ) > 

2.1 Tree level study 

The spectrum of the Dirac-Wilson operator at tree level is known and that allowed us to test 
our numerical implementation. A tree level calculation (Fig. 1(a)) shows that at fixed 6 and for 
large Nt the reweighting factor approaches the value 1, which is expected since it is a finite volume 
effect. Conversely, the variance grows as Nl increases (Fig. 1(b)). Hence a direct estimate for large 
6 angles is not reliable and we need to employ a multi-step method in order to keep the error under 
control. To this end, following [6], we factorize the relevant matrix in the following telescopic way 


g-vnAt] \ 


N-l 


Dwie)D„\0) = Y[Ae, with A^~ 1 + 0(50^) , 
1=0 



3 











Reweighting twisted boundary conditions 


A. Bussone 


0.8 


0.4 - 


- 1 -1- 1- 1- 1- 1 -1-1~ 

le-l-08 

^^^^^^^-n 

* * * 

le+07 






le-l-06 

- 





100000 

- 





3 10000 



^ 1000 



b 



100 

_ 





10 









1 


1 

0.1 

1 


0 4 8 12 16 20 24 28 32 

Nl 

(a) Mean of the reweighting factor at tree level 
for 6 = 0.1 aj a function ofNi. 
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(b) Variance of the reweighting factor at tree level 
for 6 = 0.1 oi a function ofNi. 


Figure 1: Estimates of mean and variance of the reweighting factor employing the exact formulae for the 
tree level case. Each point corresponds to a cubic lattice of the form L^. 


with Ai now near to the identity matrix as it corresponds to a small shift in the periodicity 
angle. At this point the inverse determinant and its variance are given in terms of the N analogous 
quantities, one for each multiple of 56^ in 0 , as 

1 / exp(-T7ffl’U^T7('')) \ 

detA“|J\ p(t 7 W) /p(^M)’ 


3. Monte Carlo studies 

We have calculated a number of reweighted observables in the SU(2) gauge theory with 
fermions in the fundamental representation. This theory is known to exhibit confinement and chiral 
symmetry breaking and it is therefore QCD-like. We have employed unimproved Wilson fermions 
and the Wilson plaquette gauge action. In the table below we list the configurations used in this 
work (ntcr is estimated to be —0.75(1) at j 8 = 2.2). 


V 


m 

A^cnf 

Iraj. sep. 

16 xS^ 

2.2 

-0.6 

10^ 

10 

32 X 243 

2.2 

-0.65 

394 

20 

32 X 243 

2.2 

-0.72 

380 

10 


We have reweighted observables to constant non-vanishing spatial d. To evaluate the reweighting 
factor we have used the 75 version of the Dirac-Wilson operator, since it is hermitian and with 2 
flavours if aufomafically safisfies fhe applicabilify condifion in eq. (2.1). We found fhaf for Nri > 
200 fhe reweighfing facfor on each configurafion is larger fhan 10 limes fhe average value on a few 
configuralions only. This guarantees lhal fhe average values are nol dominated by spikes and hence 
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that there are no large statistical fluctuations of the determinant. In the figures below we show the 
mean of the reweighting factor (all points are averaged over the total number of configurations). 
In Figs. 2(a) and 2(b) a good scaling for the error is clearly visible, according to and up 

to Nrf ~ 500, when increasing the number of steps. We conclude that at these parameters taking 
500—600 gaussian vectors is enough to saturate the statistical noise at the level of the gauge noise. 





(a) Mean of the reweighting factor for 6 = 0.3 as a 
function ofN-q. 





(b) Mean of the reweighting factor for 6 = ;r/2 as a 
function ofN^. 


Figure 2: Monte Carlo history of the reweighting factor, averaged over the entire number of configurations. 
The figures correspond to a volume V = 16 x 8^. 


3.1 Small volume 

Reweighting may give a significant effect in small volumes. We have first looked at the pla- 
quette because it does not explicitly depend on the BCs since, in general, in Wilson loops for each 
link entering the loop there is another one in the opposite direction and the exponential factors 
in eq. (1.1)) cancel among the two. We are neglecting autocorrelations, because measurements 
are separated by 10 to 20 molecular dynamics units. A binning procedure, including bins from 
length one up to 10, shows no significant autocorrelation effects. In Fig. 3(a) we show results after 
reweighting only one flavor, by taking the square root of the reweighting factor. No effects are 
visible in this case within errors. In Fig. 3(b) instead we have included the reweighting factor for 
both flavors and one can see a permit effect, recognizable because the plaquette is determined very 
accurately. 

The second quantity that we have studied is the pion dispersion relation. In the pion correlator 
we twisted only one flavor in the valence, hence the lowest energy state becomes a pion with 
momentum p = ±0/L. We have reweighted the correlators in order to take into account the twist 
in the sea and we have extracted the effective mass after symmetrizing the correlators in time. 

In Fig. 4 we show results for the dispersion relation and we compare them to un-reweighted 
data (i.e., with twisting only in the valence), to the continuum prediction (aS)^ = {anijff' + 30^/A^ 
and to the lattice free (boson) theory prediction cosh(a£') = 3+cosh(c!nj;i:) + 3cos(0/Az,). Over the 
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(a) Plaquette with only one flavor reweighted. (b) Plaquette with both flavors reweighted. 

Figure 3: Monte Carlo history of the reweighted plaquette. The figures correspond to a volume V = 16 x 8^. 


entire range of d values explored there is no clear effect within errors. However there is a systematic 
effect upward. The discrepancy between the reweighted data and the lattice free prediction is 
conceivably due to cutoff effects and non-perturbative dynamics. 



OINl 

Figure 4: Pion dispersion relation for V = 16 x 8^. Each point corresponds to the best determination of the 
reweighting factor. 


3.2 Large volumes 

Obviously, in large volume we do not expect the situation to improve with respect to small 
volumes, concerning both statistical accuracy and significance of the reweighting. We have looked 
at the pion dispersion relation for two rather different values of ruji. We expect to find a larger effect 
for the lighter pion but already at these volumes it appears as one can neglect the effects coming 
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from the breaking of unitarity. In Fig. 5 we show the results for the two pion masses. No sizeable 
effect is detectable in either case. 



S/Nl 

(a) Case of m ~ —0.65 which corresponds to a heavy 
pion. 



B/Nl 

(b) Case of m ~ —0.72 which corresponds to a lighter 
pion. 


Figure 5: Pion dispersion relation for V = 32 x 24^. The reweighting factor at a given 9 is obtained by a 
telescopic product involving all the previous ones. 


4. Conclusions 

We have presented an application of the reweighting method to the case of the spatial peri¬ 
odicity of fermionic fields. We have studied the approach at tree level and by mean of numerical 
simulations. For the latter we have looked at the plaquette and at the pion dispersion relation in 
small and large volumes. In both cases we have found the effects to be at the sub-percent level for 
values of 6 up to 7r/2. 
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